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CO ■ Abstract 

o 

We have calculated the Ml strength distributions in the 36 ' 38 Ar isotopes within 
large-scale shell model studies which consider valence nucleons in the sd and pf 
shells. While the Ml strength in 36 Ar is well reproduced within the sd shell, the 
experimentally observed strong fragmentation of the Ml strength in 38 Ar requires 
configuration mixing between the sd and the pf shells adding to our understanding 
of correlations across the N=20 shell gap. 
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1 Introduction 



The shell model concepts of closures, magic numbers and inert cores leads to a 
reasonably good description of nuclear phenomena in light as well as in heavy 
nuclei. However, recent experimental findings show that shell closures might 
get eroded in very neutron-rich nuclei and that the inter-shell correlations 
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Fig. 1. Calculated and experimental spectra for 36 Ar (left) and 38 Ar (right). Exci- 
tation energy E x is plotted along the y-axis. The calculations have been performed 
in the sd shell using the USD interaction [3]. The model space allows only the 
formation of positive parity states. 

among nucleons play a decisive role for the structure of these nuclei. A famous 
example is the nucleus 32 Mg which is strongly deformed in the ground state, 
despite its neutron number N = 20, e.g. [1]. Thus, shell closures are not 'rigid' 
and there have been for a long time indications for cross-shell correlations even 
in nuclei close to classical double-magic nuclei like 40 Ca. For this nucleus, these 
include the observation of a strong Ml transition [2] and of a non-vanishing 
Gamow- Teller (GT) strength built on the 40 Ca ground state [4,5,6]. On the 
theory side, recent studies [7,8] show that there are very large admixtures of sd 
configurations in low-energy states of Ca isotopes with A > 40, although these 
nuclei are usually considered as "good" pf shell nuclei. Another indication for 
the erosion of the sd-pf closure near A = 40 are the data for the Ml strength 
in 36 Ar and 38 Ar [9]. While experimentally the Ml strength distribution for 
both nuclei is rather fragmented, this feature is reproduced by sd shell model 
calculations for 36 Ar, while these studies predict only one strong Ml transition 
for 38 Ar, in clear contradiction to the data [9]. Obviously a single s<i-shell model 
space, allowing only for two proton holes with respect to a 40 Ca core, besides 
a closed neutron shell, is too limited for a description of the Ml data in 38 Ar. 

The fact that cross-shell excitations are quite important is revealed by a de- 
tailed comparison of the experimental and calculated spectrum at low energies. 
The shell model calculations have been performed within the sd shell using 
the standard USD interaction [3]. The model space does not allow the for- 
mation of negative parity states, but the calculation also misses a few of the 
low-lying + and 2 + states in 38 Ar, while the agreement between theory and 
experiment is quite sufficient for the positive parity states in 36 Ar (see Fig.l). 
In Fig. 2 we compare the experimental Ml strength distribution with the shell 
model results, which are similar to those published earlier in [9]. The shell 
model code ANTOINE [10] has been used for all calculations presented in this 
paper. 
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Fig. 2. Calculated and experimental Ml 0^~ — > lj strength distributions for 36 Ar 
(left) and 38 Ar (right). The calculations have been performed in the s<i-shell using 
the USD interaction [3]. The spin part of the Ml operator has been quenched by a 
factor of 0.75 which is in agreement with the one deduced for other nuclei in this 
mass range [11]. The same quenching of the spin part of the Ml operator has been 
used throughout this paper. Experimental error bars, which are of order of 10% for 
all measured states, are not shown here and in all further figures. 



The detailed comparison of the s<i-shell model results with the data shows 
that the Ml data for the two argon isotopes 36 Ar and 38 Ar present a challenge 
to theory: correlations between sd and pf shells should not affect the major 
features of the Ml strength distribution for 36 Ar, but must be essential for the 
38 Ar Ml data, e.g. if only two neutrons are added to the nucleus. 



The present paper is an attempt to reproduce this striking behavior on the ba- 
sis of large-scale shell model calculations which combine both sd and pf shells. 
We have followed two routes. At first, we have performed studies within the 
complete sdpf shell. However, due to computational limitations these calcu- 
lations had to be restricted to 2-particle 2-hole (2p2h) and 4p4h excitations 
from the sd- to the jo/-shell. These studies were supplemented by full diag- 
onalization calculations within the c?3/2 s i/2,/V/2P3/2 model space (denoted as 
sdpfo space below). Further additional lplh coupling between the d 5 / 2 , /s/2, 
P1/2 orbitals and the sdpfo space has been taken into account - this extended 
space is refered to as sdpfi below. As we will show in the following sections, 
the calculations in the sdpfi model space indeed result in fragmentation of 
the Ml strength for both argon isotopes and also reproduce their low-energy 
excitation spectrum. 
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2 Truncated shell model studies in the complete sdpf shells 

There have been several attempts to derive a suitable effective interaction for 
the sdpf model space (see, for example [12,13,14] and references therein). For 
our study we adopt the one of Ref. [14]. This interaction contains three parts: 
the USD interaction for the scZ-shell [3], the KB' matrix elements for the pf 
shell [15] and the G-matrix of Kahana, Lee and Scott for the cross-shell ma- 
trix elements [16]. The f 7 /2d 3 / 2 , /7/2S1/2, P3/2S1/2 andp 3 / 2 d 3 /2 monopole terms 
have been adjusted to some recent experimental data [14]. We note hat there 
is a problem with spurious center-of-mass (COM) excitations for the chosen 
configuration space. To treat this problem we have added the tenfold of the 
center-of-mass Hamiltonian to the original Hamiltonian [14]. This procedure 
pushes states with predominantly center-of-mass excitations up in energy and 
practically eliminates the spurious excitations in the low-energy states of in- 
terest here. 

Unfortunately full diagonalization studies are yet not possible in the complete 
sdpf model space due to computational limitations. However, it is feasible 
to perform calculations for 38 Ar in which 2 particles or even 4 particles are 
promoted from the sd shell to the pf shell. We will refer to these studies as 
2p2h and 4p4h calculations, respectively. The dimension of the 4p4h model 
space is 227 625 436 for M = + states in the m-coupling scheme. For 36 Ar we 
have performed 2p2h calculations, corresponding to 3 148 356 M = + states. 

These studies, however, showed no improvement for the 38 Ar spectrum com- 
pared to the pure scZ-shell results. In fact the first excited + state lies at an 
excitation energy above 7 MeV, with very little mixing to the ground state 
which had basically an sd shell configuration. Consequently no fragmentation 
has been found in our calculation of the Ml strength distributions for 38 Ar. 
We also have performed No-Core Shell Model [17,18] calculations in a 4hu 
model space using the UCOM interaction [19] obtaining similar results. 

The mixing of the sd and pf shell configurations can be enhanced by slight 
monopole modifications in the interaction. This strategy has been applied in 
large-scale shell model calculations of 36 Ar and leads to a quite satisfying de- 
scription of the strongly deformed 4p4h band and its coupling to the ground 
state band [20]. Following this reference we have lowered the energy of the 
npnh configurations with n = 4 and n = 6 by suitable monopole shifts, in 
this way enforcing p/-shell configurations to mix into the ground states. As 
it is expected the mixing results in fragmentation of the Ml strength distri- 
bution. This can be seen in Fig. 3 which shows the calculated Ml strength 
for 38 Ar obtained in a 4p4h calculation and employing Lanczos diagonal- 
ization with 50 iterations. Indeed the strong Ml transition observed in the 
pure s<i-shell calculation is now split into several components and some Ml 
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Fig. 3. Calculated (upper panel) and experimental (lower panel) Ml strength dis- 
tributions in 38 Ar in the full sdpf space with 4p4h cross-shell excitations. The cal- 
culations have been performed using a monopole-modified interaction as described 
in the text. 

strength is shifted to higher excitation energies. Nevertheless, the agreement 
with data is only marginal. In particular, the data show already some Ml 
strength around 8 MeV, which is missing in the calculation. This is possi- 
bly improved if npnh configurations with n > 4 were included - previous 
shell model studies [7,8,13,14] indicated that energies of the states in nuclei 
around 40 Ca are converging rather slowly when npnh cross-shell excitations 
with larger n are gradually added. While such converged calculations are yet 
impossible in the complete sdpf shell, they are, however, feasible in the sdpf 
and sdpf i model spaces (see text above) to which we turn in the next section. 



3 Large-scale shell-model calculations in the truncated sdpf model 
space 



We have performed diagonalization calculations in the sdpfi model space 
which includes the spin-orbit partners (missing in the sdpfo space) approx- 
imately by allowing lplh excitations from the c/5/2 orbital to the rest of the 
sd shell, and from the /7/2P3/2 space to the /5/2P1/2 sub-space, respectively. 
In other words we put lplh predominantly spin-flip excitations on top of the 
complex npnh (si/ 2 d 3 /2) A ~ 28 ~~ n (f7/2P3/2) n cross-shell configurations obtained 
within the sdpfo space. This extension increases the dimension of the model 
space from 2.3 million (sdpfo) to nearly 75 million (sdpfi) for 38 Ar. For the 
calculations we have used the same effective interaction as in [7]. The sdpfo 
model space allows up to 8 or 10 particles to be promoted from the s<i-shell 
to the pj-shell for 36 Ar and 38 Ar, respectively. It turns out these cross-shell 
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Fig. 4. Calculated and experimental spectra for 38 Ar within the full sdpfi space. 



excitations are sufficient to achieve a good agreement between our calcula- 
tions and the experimental 38 Ar spectrum. As one can observe in Fig. 4, the 
calculation now yields the low-lying + and 2 + excitations which were miss- 
ing in the sd-model space (see Fig. 1). The calculated 36 Ar spectrum in the 
sdpfi model space (see Fig. 5) agrees similarly well with the data as the pure 
sa?-shell result (Fig. 1) for the positive parity states. For both argon isotopes 
the calculated negative parity spectrum, missing in the pure sd-shell studies, 
agrees quite well with the observed one. 

It is interesting to inspect more closely the role of npnh cross-shell and lplh 
spin-flip excitations in the sdpfi space for the 36 Ar and 38 Ar ground states. 
To quantify the importance of npnh cross-shell excitations we have calculated 
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Fig. 5. Calculated and experimental spectra for 36 Ar within the full sdpfi space. 
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Table 1 

The structure of the ground state and of two final states with strong Ml transitions 
(the 1 + , T = 1 at 10.14 MeV and the 1+, T = 1 at 14.77 MeV states) for 36 Ar in the 
full sdpfi space. The coefficient A n (n d f) gives the contribution of the component 
with n particles promoted across the sd-pf shell gap. The n d f quantity indicates 
whether lplh spin-flip components are present (n d f = 1) or not (n d f = 0). Only- 
components with weights larger than 0.01 are shown. 
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the weights A n (n d f = n d + rif) of the 

(4/2) 6 ^(s 1 / 2 rf3/2) A ~ 28 - ri+rid (/7/2P3/2) ri -^(/ 5 /2Pl/2)^ 

configurations, where n d f = n d + n,f define the total number of excitations from 
the d 5 / 2 orbital to the rest of the s<i-shell and from the (/7/2P3/2) subspace to 
the respective spin-orbit partners. Our calculations are restricted to n d f = 1. 
Thus, with A n (n d f = 1) we have identified the components of the wave func- 
tions which arise from lplh couplings to the spin-orbit partners within the 
sd-shell (rid = 1,?V = 0) or p/-shell {n d = 0,n/ = 1), respectively. The dom- 
inant configurations of the ground state wave functions are given in Tables 1 
and 2 for 36 Ar and 38 Ar, respectively. In the sdpfi model space, approximately 
half of the 36 Ar ground state corresponds to pure s<i-shell configurations, where 
the lplh excitations of the d$/2 orbit play a tiny role (1%). 

Cross-shell 2p2h components contribute about 41% to the 36 Ar ground state, 
where the couplings to the spin-orbit partners are now sizable (11%). Higher- 
order cross-shell components with n > 2 are relatively unimportant (less than 
10%). These components, however, are larger in the 38 Ar ground state. We 
find that this state is dominated by npnh cross-shell excitations with n < 4, 
with OpOh and 2p2h excitations amounting to about 75% of the wave function. 
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Table 2 

The structure of the ground state (0 + ,T = 1) and of two final states with sizable 
Ml transitions (the 1+, T = 1 at 12.96 MeV and the 1+, T = 2 at 18.22 MeV states) 
for 38 Ar in the full sdpfi space. The coefficient A n {n<if) gives the contribution of the 
component with n particles promoted across the sd-pf shell gap. The ridf quantity 
indicates whether lplh spin-flip components are present (ridf = 1) or not (ridf = 0). 
Only components with weights larger than 0.01 are shown. 
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Importantly, the pure sd-shell configuration contributes only 33% to the 38 Ar 
ground state. The rest corresponds to excitations to the pf shell and hence 
opens new channels for the fragmentation of the Ml strength distribution, to 
which we turn now. 

We have calculated the Ml strength distributions by using the Lanczos method 
with 200 iterations. This is sufficient to achieve convergence for the calculated 
energies of the lowest states. However, at higher excitation energies the cal- 
culated Ml distributions represent only the total strength per energy interval 
rather than true states. We note that we have scaled the spin part of the Ml 
operator by a constant factor q = 0.75, as it is customary in shell model calcu- 
lations. It has been shown that shell model calculations using such an effective 
spin operator describe the Ml [11] and GT [21,22] transitions of ^/-nuclei quite 
well. 

The calculated Ml strength distributions for 36 Ar and 38 Ar within the full 
sdpfi model space are shown in Fig. 6. For 36 Ar the result is similar to the 
srf-shell calculation as, in agreement with the experimental data, it shows a 
rather strong transition at around 10 MeV. At higher energies, the present 
36 Ar Ml strength distribution is significantly more fragmented than the one 
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Fig. 6. Comparison of the experimental and theoretical Ml strength distributions 
for 36 Ar (left) and 38 Ar (right) calculated in the sdpfi model space. 

obtained within the pure s<i-shell. Our calculation also predicts a noticeable 
amount of Ml strength outside of the experimental energy window. In fact, we 
calculate a total Ml strength for 36 Ar of 6.76 fi 2 N . Experimentally a summed 
Ml strength of 2.65(12) /j 2 N is observed up to an excitation energy of 13.8 
MeV, which is in agreement with the calculated value (2.8 ^ 2 N ). However, our 
calculation predicts an Ml strength of 1.4 ijl 2 n in the energy interval 13.8-15 
MeV, where experimentally no strength has been observed. 

For 38 Ar the result from the sdpfi model space is strikingly different than that 
of the pure s<i-shell calculation reported above. The calculated Ml strength 
in the energy regime below 15 MeV is now strongly fragmented like the ex- 
perimental one. Moreover, we find a total Ml strength of 2.88 jj, 2 N in the 
experimental energy window, to be compared with the experimental value of 
2.86(18) fi 2 N . For 38 Ar (with a T = 1 ground state) Ml transition can lead 
to final states with T = 1 and T = 2. Most of the Ml strength to T = 2 
states resides at excitation energies above 15 MeV so that the main peak of 
the Ml strength distribution around E x = 13 MeV corresponds dominantly 
to AT = transitions. As for 36 Ar, our calculation predicts a strong part of 
the Ml strength to reside at energies higher than the current experimental 
energy window. In fact, our calculated total Ml strength is 6.01 ji 2 N for 38 Ar. 
It would be desirable to experimentally search for this predicted Ml strength 
at higher energies. 

To explore the origin of the fragmentation observed in the Ml distributions for 
both argon isotopes, we have also analyzed the wave functions of the strongest 
1 + states. 

For 36 Ar, Table 1 lists the dominant configurations of the 1 + ,T = 1 state at 
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Fig. 7. Evolution of the Ml strength distribution with increase of the model space. 
Here nt = n + ridf , where n is the number of particles excited across the sd-pf gap 
and n<if < 1 stands for the number of nucleons promoted from the d§/2 orbital to 
the si/2^3/2 subspace and from the /7/2P3/2 to the /5/2P1/2 ones. 

10.14 MeV, which corresponds to the single strong Ml transition observed in 
the calculation (and in the data), and the 1 + , T = 1 state at 14.77 MeV, which 
is a typical state in the energy range of strongly fragmented Ml transitions. 
The 1 + state at 10.14 MeV is dominated by OpOh and 2p2h configurations 
(> 80%), similarly to the ground state. We notice that the state has a quite 
significant 2p2h component (40%) with a simultaneous cross-shell excitation 
of a proton and a neutron. Thus, we find, as in the pure s<i-shell calculations, 
a strong Ml transition around 10 MeV, however, the configurations of the 
initial and final states are quite different in the sd and the sdpf\ model spaces. 
The 1 + state at 14.77 MeV has sizable 4p4h configurations (about 30%), but 
in contrast it has a rather small OpOh component (less than 10%). In this 
energy range the various configurations are quite strongly mixed resulting in 
the strong fragmentation of the Ml strength. 



For 38 Ar Table 2 lists the dominant configurations of the 1 + ,T = 1 state at 
12.96 MeV which carries a sizable fraction of the fragmented Ml strength 
below 15 MeV. Different to the ground state, OpOh excitations (pure sd-shell 
configurations) amount only to 5%, while 2p2h and 4p4h excitations are dom- 
inant. But importantly we also find a noticeable contribution of 6p6h configu- 
rations in these states. It is further worth mentioning that neutron excitations 
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are more important than proton excitations. 



Obviously there are many 1 + states with large configuration mixing at mod- 
erate excitation energies which have a substantial p/-shell contribution, but 
only a rather small pure sd-shell component or a nearly vanishing OpOh com- 
ponent. Table 2 lists the dominant configurations of the 1 + ,T = 2 state at 
18.22 MeV, which is a typical example in this energy regime. 



Finally we have explored the evolution of the Ml strength distribution in 
38 Ar with the number of allowed cross-shell excitations. It is convenient to 
characterize the different truncations using the quantity n t = n + n d f, where 
n is the number of particles moved across the sd-pf shell gap and n d f in- 
dicates the number of nucleons moved from the d 5 / 2 orbital to the Si/ 2 d 3 /2 
subspace and from the /7/2P3/2 subspace to the spin-orbit partners /5/2Pi/2- 
In our sdpfi calculations we have allowed only lplh excitations between spin- 
orbit partners, hence n d f < 1. The simplest excitation scheme is n t = 2 which 
allows the mixing of three configurations, namely the one with n = 0, n d f = 0, 
n = 2,n d f = and n = 0,n d j = 1. (Note that n = 1 excitations are not 
possible due to parity.) The n t = 2 excitations do not change qualitatively 
the shape of the Ml strength function (see Fig. 7) as compared to the pure 
sa?-shell results. Only the centroid is moved to higher energies, as the origi- 
nal single spin-flip state is shifted up in energy and is connected by the Ml 
operator only to the n = 0, n d f = configuration. If the number of excita- 
tions is increased to n t = 3, only one new configuration with n = 2, n d f = 1 
is added. This configuration is, however, crucial for the fragmentation of the 
Ml strength (see Fig. 7). Besides the large Ml matrix element between the 
n = 0, n d f = and n = 0, n d f = 1 configurations in the s<i-shell, there is now 
another big Ml transition in the model space (between the n = 2, n d f = and 
n = 2, n d f = 1 configurations) involving 2p2h excitations across the shell gap. 
The interference of these two contributions produces the fragmentation of the 
Ml strength. The further addition of new components like n = 4, n d f = 
(n t = 4) and n = 4, n d f — 1 (nt — 5) does not change the fragmentation 
picture qualitatively, however, the summed B(M1) strength increases from 
2.95 fi 2 N (n t = 3 case) to 4.02 n 2 N and 4.20 n 2 N , respectively, as more p/-shell 
components are mixed into the ground state wave function. We also observe 
that the convergence for the energies of many states is not even achieved at 
the nt — 5 level (see Fig. 8). This indicates that more complicated configura- 
tions with n t > 5 have to be included, which we have done in our complete 
n = 10, n d f — 1 (n t — 11) calculation discussed above. 
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Fig. 8. Comparison of experimental and calculated spectra for Ar in the sdpfi 
space for different values of rit = n + ridf with ridf fixed at 1 and varying n cor- 
responding to npnh cross-shell excitations. The full sdpfi space calculation corre- 
sponds to n t = 11. 

4 Conclusion 



We have performed large-scale shell model calculations for 36 Ar and 38 Ar 
within the complete 51/2^3/2/7/2^3/2 model space and additionally allowing 
lplh excitations from d 5 / 2 to the si/ 2 , d 3 / 2 orbitals or from f 7 / 2 and p 3 / 2 to 
the rest of the p/-shell. We reproduce the energy spectra for both isotopes, in- 
cluding the negative parity states. Importantly our calculations identify sizable 
contributions of cross-shell components in the 36 Ar and 38 Ar ground states, 
implying the onset of erosion of the N=20 shell closure. However, this erosion 
is a slow process involving many npnh excitations until, for example, conver- 
gence of the energies of the lowest excited states is achieved. These studies 
reveal the onset of fragmentation in the Ml strength distributions. However, 
only a marginal agreement with the data has been obtained, which might be 
improved if higher-order cross-shell excitations are considered. 
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These cross-shell excitations play an important role for the Ml strength dis- 
tributions. For 36 Ar our calculation exhibits a single strong Ml transition at 
10 MeV, which is in agreement with the measured distribution. At higher ex- 
citation energies the calculated distribution is quite strongly fragmented, in 
contrast to the results of pure sd-shell calculations. Most of this fragmented 
strength resides, however, outside of the energy range for which currently data 
exist. 

For 38 Ar our calculated Ml distribution is in even stronger contrast to the 
one obtained in the pure srf-shell. The cross-shell npnh excitations lead to 
a very strong fragmentation of the strength. Particularly important are here 
( s i/2^3/2) 10 ~ n (/7/2P3/2) n configurations with n = 2,4 and 6 which are strongly 
mixed in the wave functions of 1 + states above 10 MeV. Consequently, also the 
configurations which include lplh spin-flip excitations involving the d 5 / 2 , fs/2 
and pi/2 partners build on top of the npnh cross-shell excitations get strongly 
mixed as well. As a consequence the Ml strength distribution is strongly 
fragmented, as observed in the data. As for 36 Ar our calculation predicts also 
for 38 Ar a sizable amount of Ml strength to reside at energies higher than the 
current observational limit. 

Very recent experimental data [24] obtained via the 40 Ar(7, 7') photon scat- 
tering reaction suggest a very strong fragmentation of the Ml strength in 
40 Ar, where only one 1 + state with a B(M1) value of 0.145(59) fi 2 N has been 
identified in the energy region between 7.7 and 11 MeV. The evolution of the 
Ml strength with increase of the neutron number (N > 20) for Ar isotopes is 
an intriguing issue which will contribute further to our understanding of the 
N = 20 cross-shell dynamics. 
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